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Abstract 

Background: Ubiquitin-mediated protein modification and degradation are believed to play important roles in 
mammalian spermatogenesis. The catalogues of ubiquitin activating enzymes, conjugating enzymes, and ligases 
(E3s) have been known for mammals such as mice and humans. However, a systematic characterization of E3s 
expressed during spermatogenesis has not been carried out. 

Results: In present study, we set out to mine E3s from the mouse genome and to characterize their expression 
pattern, subcellular localization, and enzymatic activities based on microarray data and biochemical assays. We 
identified 398 putative E3s belonging to the RING, U-box, and HECT subfamilies and found that most genes were 
conserved between mice and humans. We discovered that 73 of them were highly or specifically expressed in the 
testes based on the microarray expression data. We selected 10 putative E3 genes to examine their mRNA 
expression pattern, and several genes to study their subcellular localization and E3 ligase activity. RT-PCR results 
showed that all the selected genes were predominately expressed in the testis. Some putative E3s were localized in 
the cytoplasm while others were in both the cytoplasm and the nucleus. Moreover, all the selected proteins were 
enzymatically active as demonstrated by in vitro and in vivo assays. 

Conclusions: We have identified a large number of putative E3s that are expressed during mouse spermatogenesis. 
Among these, a significant portion is highly or specifically expressed in the testis. Subcellular localization and 
enzymatic activity assays suggested that these E3s might execute diverse functions in mammalian spermatogenesis. 
Our results may serve as an initial guide to the field for further functional analysis. 
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Background 

Spermatogenesis, a complex yet highly regulated develop- 
mental process of male germ cells, consists of three 
stages— the proliferation of spermatogonia, the meiosis of 
spermatocytes and the morphogenesis of spermatids [1,2]. 
The formation of the terminally-differentiated germ cells, 
the spermatozoa, requires a series of important processes 
that are unique to spermatogenic cells, including nuclear 
condensation, mitochondrial rearrangement, histone re- 
placement by transition proteins and protamin, the shed- 
ding of residual bodies, and the formation of acrosomes 
[3]. It is not surprising that a large number of testis- 
specific genes are required for these events. In mice and 
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rats, several high throughput gene expression datasets 
have revealed that a large number of genes are expressed 
to fuel spermatogenesis [4,5]. It is also estimated that 
about 4% of protein-coding genes in the mouse genome 
are specifically expressed in the mouse testis [6]. 

Ubiquitination of protein is an indispensable post- 
translational modification that serves as a component of 
the protein quality control system. It is also involved in 
diverse biological processes such as signal transduction, 
DNA repair, transcriptional regulation in a protein 
degradation-independent way [7]. In principle, ubiquiti- 
nation is a process containing three ubiquitin transfer- 
ring reactions catalyzed by three corresponding 
enzymes. At first, the 8.5-kDa ubiquitin polypeptide is 
activated by the ubiquitin- activating enzyme (El) in an 
ATP-dependent way. Subsequently, this activated ubiqui- 
tin is transferred to the ubiquitin conjugating enzyme 
(E2) via a thioester bond. In the end, the ubiquitin is 
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transferred to a lysine residue of the substrate catalyzed 
by the ubiquitin-protein ligase (E3) [8]. Substrate specifi- 
city is conferred by E3s, implying a much larger number 
of E3s in the genome than Els and E2s. The mammalian 
genome encodes 1-2 Els, 10 ~ 20 E2s and several hun- 
dred E3s [9]. There are three typical E3 families includ- 
ing the RING finger family, HECT domain family and 
the U-box family [10-12]. 

To our knowledge, approximate 29 E3s have been 
reported to be expressed in mammalian testes and most 
of them execute diverse functions at different stages of 
spermatogenesis (summarized in Additional file 1). A 
dozen of E3s as well as associated complex play import- 
ant roles in DNA double strands break (DSB) and his- 
tone modifications during early meiosis. As examples, 
the expression of LASU1 complex and Cul4A-CRL4 is 
detected primarily in spermatogonia stage and involved 
in histone ubiquitination during meiosis [13-16]. An N- 
end rule pathway E3 ligase, UBR2, plays a critical role in 
transcriptional chromosome inactivation via ubiquitina- 
tion of histone H2A [17]. The majority of E3s that are 
highly expressed in haploid germ cells probably contrib- 
ute to the morphogenesis during the formation of 
spermatozoa. For instance, MARCH 10 is a microtubule 
associated E3 and is involved in the organization and 
maintenance of the flagella of spermatozoa [18]. 

Although it has been known that the mammalian gen- 
ome encodes a large number of E3s, the exact number is 
still not known particularly when higher specificity of 
mining is a major concern to the experimental scientists. 
Moreover, how these putative E3s are expressed during 
mammalian spermatogenesis is still an open question. In 
the present study, we mined out putative E3s from the 
mouse and human genome, evaluated their expression 
in multiple tissues, particularly in the mouse testis at dif- 
ferent stages. The ligase activity of selected E3s was con- 
firmed by in vitro and in vivo assays. Our list of E3s 
expressed during spermatogenesis provides a valuable 
source for future functional studies of the ubiquitination 
during mammalian spermatogenesis. 

Results 

Mining of E3s from the mouse and human genomes 

To mine putative E3s from the mouse and human gen- 
omes, we first compiled all protein-coding genes from sev- 
eral microarray datasets (GSE96, GSE97, GSE1133, 
GSE2361,GSE9954, GSE1986) generated from gene ex- 
pression profiling of multiple tissues and the EST dataset 
from the UniGene database [19-22]. As a result, 26762 
and 23058 mouse and human genes were identified, re- 
spectively, and 15952 genes were homologous genes 
(47.1%) based on the HomoloGene database annotation 
(Figure 1A). We then searched all protein coding genes 
for domains/motifs in the Pfam_ls and Pfam_fs library 
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Figure 1 Numbers of putative E3s mined from the mouse and 
human genomes and the homologous genes between the two 
species. Gene list of protein-coding genes was compiled from lists 
of genes that can be detected using the Affymetrix microarray 
technology. Sequences of these genes were downloaded from NCBI 
Refseq database and were scanned for RING/U-box/HECT domains 
using the HMMER software. Homologue identification was based on 
the annotation of NCBI HomoloGene database. (A) Number of 
homologues of all protein-coding genes in humans and mice. The 
number in the intersection is the number of homologues while the 
other two represent genes unique to the two species. (B) Number 
of homologues of putative E3 genes in humans and mice. (C) Numbers 
of different E3 genes in humans, mice and yeasts. 



(release 23.0) using HMMER 2.3.2 software package 
(http://hmmer.janelia.org/) with most of the parameters 
being set to the default values. The E value of the hit was 
set to be no more than 0.1. Proteins containing the RING 
domain (zf-C3HC4), the HECT domain and the U-box 
domain were considered as typical E3s. As a result, 398 
and 411 putative E3s were identified from the mouse and 
the human genome, respectively (see Additional file 2 for 
mouse E3s). Among them, 335 putative E3s were homolo- 
gues (70.7%) between the two species (Figure IB, Add- 
itional file 3). The proportion of homologues in the E3 set 
is significantly higher than expected from the number of 
general homologues of the mouse and the human gen- 
omes (P-value is 3.8E-29 based on the binomial test). We 
also identified the yeast homologues of the initial 335 
mouse/human homologous pairs from 50 yeast E3s 
(Figure 1C), and only 7 highly conserved E3s among these 
three species were found (Table 1). These observations 
suggested that E3 genes are more conserved between the 
mouse and human genome. 

Expression of E3 mRNAs in the mouse testis 

We next examined the mRNA expression levels of E3 
genes in the mouse based on four microarray datasets 
(GSE97, GSE1133, GSE9954, and GSE1986) and the Uni- 
Gene dataset. For each dataset, we defined 5 expression 
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Table 1 E3s conserved in the yeast, the mouse and the human 
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levels of an mRNA in a particular tissue in the increasing 
order of tissue specificity— A (absent), P (present), HP 
(highly present), MS (multiple tissue-specific), and SP 
(specific). The level A/P is just based on a genes P/A call 
value determined by the Affymatrix platform using the 
MAS5 algorithm [23]. The other three levels were based 
on the comparison of the mRNAs z-score with three 
threshold values (See Methods for further explanation). 
The evaluations from different datasets were summarized 
by voting— A level value is assigned only when it is sup- 
ported by no less than two datasets. A putative E3 is al- 
ways assigned to the highest specificity group exclusively. 
As a result, we identified 267 P level E3s, 39 HP level E3s, 
2 MS level E3s, and 32 SP level E3s in the mouse testis 
(Table 2). Compared with other tissues such as the liver, 
the lung, the muscle, the testis contains significantly more 
specific E3s (Figure 2). Therefore, a large number of E3s 
are expressed in the mouse testis and a significant portion 
is specifically expressed. 

Stage-specific expression of E3 mRNAs during mouse 
spermatogenesis 

We further studied the stage-specific expression of puta- 
tive E3s during mouse spermatogenesis using three micro- 
array datasets (GSE640, GSE926 and GSE4193) derived 
from gene expression profiling of testis at different days 
postpartum or of different spermatogenic cell types 
[6,24,25]. Genes that are 3-folds or more differentially 
expressed at one or more time point (or cell type) were 
selected for further clustering analysis. Based on visual in- 
spection of the heat maps of the clustering analysis, three 
clusters of genes were identified according to the stages 



Table 2 Numbers and mRNA expression levels of 
annotated E3s in the mouse testis 
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when their expressions start to increase. Genes in group 
MI are expressed at a higher level at the mitotic division 
stage, genes in group ME start to increase their expression 
level at meiotic division stage while genes in group PM 
start to increase their expression at the post meiosis stage. 
Results from the three datasets were again summarized by 
voting— a gene was assigned to a particular group only if it 
was supported by 2 or 3 votes. As results, 20, 24 and 6 E3s 
were assigned to the MI, ME, and PM groups, respect- 
ively, totaling up to 50 (Table 3 and Additional file 4). By 
taking intersections of each of three groups with each of 
the SP, MS, HP, and P classes, it can be seen that the puta- 
tive testis specific E3s increase their expression during or 
after meiosis. 

Stage-specific expression of E3 proteins during mouse 
spermatogenesis 

In a separate study, we examined protein levels in type 
A spermatogonia, pachytene spermatocytes, round sper- 
matids, and elongated spermatids using the iTRAQ mass 
spectrometry (unpublished data). Thousands of proteins 
were identified in each cell type. We selected E3s from 
these detected proteins and examined their dynamic 
changes during spermatogenesis. 38 annotated E3s were 
found (Additional file 5). The heat maps of the cluster- 
ing analysis showed that these proteins could be sepa- 
rated into two groups according to the dynamic change 
of their expression. Protein expression of group 1 is up 
regulated after meiosis while group 2 is on the contrary. 
Meanwhile, we performed an RNA deep sequencing in 
the same cell type (unpublished data). Compared to 
iTRAQ analysis, the expression patterns of most RNAs 
were consistent with those of protein expression. The 
changes of a portion of proteins were lagging behind 
those of their RNAs probably due to translation inhib- 
ition (Figure 3). 

Validation of mRNA expression of selected putative E3s 
by RT-PCR 

Among the 32 putative testis-specific E3s, 19 have not 
been well studied. We selected 10 of them randomly and 
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Figure 2 Expression profiles of E3 genes in different tissues. Expression of different families of E3s were ranked based on their levels in the 
mouse testis (P for present, HP for highly present, SP for specific). The expression of genes of these three levels were examined in other tissues 
such as the liver, the lung, and muscle. (A) Almost all the SP-level E3 genes were only expressed in the mouse testes. (B) More HP-level genes 
were expressed in the testes than in other tissues. (C) Genes that were expressed in the testis (P-level genes) were also expressed in other tissues. 





checked their mRNA expression in 10 different tissues 
(testis, ovary, uterus, liver, lung, muscle, heart, brain, 
spleen and kidney) using RT-PCR. As shown by 
Figure 4A, majority were exclusively detected in the 
testis while the remaining one were expressed in the 
testis at much higher levels than in other tissues. These 
result indicated that our testis specificity evaluation of 
genes based on microarray data is highly reliable. Exam- 
ination of their expression at different days postpartum 
indicated that almost all of them were expressed at 
stages when haploid cells are generated (Figure 4B). We 
further examined expression of some genes in isolated 
type A spermatogonia (SG-A), pachytene spermatocyte 
(pacSCs), round spermatid (rSTs) and elonged spermatid 
(eSTs). The purity of the isolated germ cells all exceeded 
90% as determined by morphological evaluation, and 
was also confirmed by measuring the expression of 
known marker genes that are either uniquely or highly 
expressed in each type cells (Gfral and Sohlhl for type 
A spermatogonia; Ldhc for germ cell except spermato- 
gonia; Sycpl, Piwil2 and Hit for spermatocytes; Tnpl 
and Prm2 for spermatids; Fshr for Sertoli cell). As shown 
by Figure 4C, most of the E3s start their expression in 
either pacSCs or rSTs. 



Table 3 Numbers of stage-specifically expressed E3 
mRNAs during mouse spermatogenesis 
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Subcellular localization of putative E3s with or without 
transmembrane domains 

A significant number of E3s have signal peptides and/or 
transmembrane domains (TMD) suggesting that they 
may anchor to the cells membrane system to execute 
their functions. 16 and 65 annotated E3s were identified 
with either signal peptide and/or TMD, respectively 
(Figure 5A). As subcellular localizations of E3s help to 
understand their function, and no antibodies were avail- 
able for in vivo studies, we decided to transfect CHO 
cells with plasmid constructs of E3s fused to EGFP to 
study their subcellular localization. We selected 2 E3s 
with TMD (RNF133, ZNRF4) and 3 without (RNF151, 
RNF126, NSMCE1), and used mouse Ubc6, an endoplas- 
mic reticulum-localized integral membrane E2, as the 
reference [26]. CHO cells were co-transfected with the 
EGFP-E3 and the MmUbc6-DsRed constructs. Subcellu- 
lar localization of the fusion proteins were inspected 
under fluorescent microscope. As shown in Figure 5B, 
RNF151 was only localized in the nucleus in a punctate 
manner. RNF126 and NSMCE1 was localized both in 
the nucleus and the cytoplasm in a diffused manner, 
while RNF133 and ZNRF4 were only localized in the 
cytoplasm, co-localized well with MmUbc6-DsRed. The 
localizations of RNF133 and ZNRF4 were consistent 
with previous report [27,28], indicating that our assays 
using the fluorescent fusion proteins were reliable for 
subcellular localization studies. 

In vitro and in vivo ubiquitin ligase activities of putative 
E3s 

It was important to test whether the putative E3s indeed 
possess ubiquitin ligase activities in order to understand 
their function. However, due to the lack of antibodies to 
immunoprecipitate native E3s from testis or germ cell 
lysate, we decided to used recombinant protein and 
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Figure 3 Putative E3s whose proteins were detected in one or more germ cell types using the iTRAQ mass spectrometry and their 
corresponding RNA levels detected by RNA-seq data. The protein and RNA expression data of these putative E3s were from two 
independent studies aiming to examine global gene expression at the mRNA and protein levels. Note most putative E3s were expressed at or 
after meiosis. The mRNA expression either matches or precedes their protein expression. SG-A: type A spermatogonia, pacSC: pachytene 
spermatocyte, rST: round spermatid, eST: elonged spermatid. 



tagged protein immunoprecipitated from transfected 
cells to examine their enzymatic activity. Three putative 
E3 genes (Nsrncel, Rnfl26 and RnflSl) were cloned into 
the pGEX-4Tl vector. GST fusion proteins were 



expressed and purified from E coli. In vitro ubiquitina- 
tion assays were performed by including El, E2, E3, and 
ubiquitin. Anti-Ubiquitin (Ub) antibody was used to de- 
tect polyubiquitin chains. To screen the E2 with highest 
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Figure 4 Examination of mRNA expression of annotated murine testis-specific E3s by RT-PCR. (A). mRNA expression in 10 tissues. 
(B). mRNA expression in testis of different days post partum. (C). mRNA expression in 4 isolated germ cell types. 
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Figure 5 Analysis of E3s containing signal peptide and transmembrane domains (TMDs). (A) Number of E3s which possessed signal 
peptides and TMDs. (B) CHO cells were transiently transfected with GFP tag recombinant plasmids. The cells were then processed for 
immunofluorescence microscopy 24 hours after transfection and MG-132 treatment. The scale bar is 50 urn. 



ubiquitin catalytic activity, three human and one mouse 
recombinant E2s (UBE2D2/human, UBE2D3/human, 
UBE2L3/human and Ube2k/mouse) were tested. Our 
results showed that UBE2L3 has barely-detectable activ- 
ity while UBE2D2 has the highest activity (Figure 6). 
The polyubiquitin chain was formed when all reactants 
were added (El, E2, E3, and ubiquitin), while, it was not 
detected if any of the reactants were excluded. There- 
fore, the ligase activity of the putative E3s was detected 
in a substrate-independent way. 

The ligase activity of these E3s was also confirmed in 
HEK 293FT cells transfected with plasmid constructs 
that gave rise to FLAG -tagged E3s and HA-tagged Ub. 
Cell lysate was subjected to immunoprecipitation (IP) 
using the FLAG antibody conjugated beads and immu- 
noblotted using the anti-HA antibody. As shown in 



Figure 7A, poly (Ub) chains were detected for all tested 
E3s. It has been known that some E3s are first polyubi- 
quitinated before they transfer the poly (Ub) chain to 
their substrates. To tell whether the poly (Ub) chain 
detected was from the E3s or any protein that was co- 
immunoprecipitated with the E3s, the cell lysate was 
first boiled with high concentration SDS to disrupt any 
potential protein-protein interaction before immunopre- 
cipitated with the beads. As shown in Figure 7B, polyu- 
biquitination was still detected with RNF126 and 
RNF151. These results showed that some E3s promoted 
self-ubiquitination before their substrates were ubiquiti- 
nated. A potential role of self-ubiquitination of E3s is to 
regulate the ligase activity and recruit substrates during 
spermatogenesis [29]. Our results also indicated that it 
was very likely that most, if not all, E3s identified by 
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Figure 6 Ligase activity of putative E3s was validated by in vitro assays. Polyubiquitin chains were formed in the presence of UBE2D2 and 
Ube2k. Polyubiquitin chains were observed only in the presence of E1, E2, E3 and ubiquitin. GST protein was used as negative control (right 
Iane6). 



bioinformatics analysis were active enzymes. Considering 
their different sub-cellular localizations as indicated by 
the transfection assays, it seemed that these E3s were 
involved in diverse molecular processes during sperm- 
atogenesis such as protein quality control, organelle 
turnover, chromatin remodeling, to name a few. 

Discussion 

In the present study, we conducted genome-wide search 
and annotation of putative E3s from several high- 
throughput expression datasets. We mined out putative 
E3s from protein coding genes whose expression have 
been examined in different tissues, and investigated their 
mRNA and protein expressions in the mouse testis. 
According to microarray data, we discovered that a large 
number of putative E3s were expressed in the mouse 
testis and that some of them were specifically expressed 
there, which were confirmed by RT-PCRs. The protein 
expression patterns of 38 E3s were detected in different 
spermatogenic cells using the iTRAQ MS technology 
(our unpublished data). The E3 ligase activity of some 
selected proteins was confirmed by ubiquitin assays. 
These results indicated that our mining and annotation 
strategy was reliable and distinguished its self from some 
previous mining procedures purely based on sequence 
analysis. 



According to a report by Li et al. based on bioinfor- 
matics analysis, the human genome contains 617 puta- 
tive E3s [9]. This number was derived from motif 
searching using the Hidden Markov Models (HMMs) for 
domains such as HECT, RING/U-box, F-box, SOCS-box 
etc. It is known that most, if not all, HECT and RING/ 
U-box proteins are true E3s while those F-box-, SOCS- 
box containing ones are not necessarily a component of 
the E3 complexes [30]. Therefore, we limited our 
searches by using only the HECT and RING/U-box 
domains to reduce the false positive rate. Indeed, Li 
et al. used several datasets containing predicted tran- 
scripts and corresponding proteins of the human gen- 
ome in order for the search to be comprehensive. 
However, we think for bench-working scientists, higher 
specificity is more important than higher sensitivity. Ac- 
cordingly, we only used protein-coding genes derived 
from the Affymetrix expression GeneChips (Mouse Gen- 
ome 430 2.0 Array, Murine Genome U74 Version 2 
Array, Human Genome U95A Array, Human Genome 
U133A Array and GNF1M/GNF1H custom array) and 
the UniGene database as the target sets to be mined. As 
a result, we identified 398 putative E3s, of which 360 
were the RING family members from the mouse gen- 
ome. Our 411 human E3s are not all covered by Lis 617 
E3s. The intersection between our 411 human E3s and 
the 617 by Li et al. is 308. We found that the list of Li 



Hou et al. BMC Genomics 2012, 13:495 
http://www.biomedcentral.eom/1 471 -21 64/1 3/495 



Page 8 of 12 



IP: FLAG 



IB: FLAG 



IB: HA 



FLAG 



170- 
130- 
95- 
72- 
55- 
43- 
34- 
26- 



118: 

95 
72- 
55 
43- 
34- 
26- 



RNF126 



118: 

95- 
72- 
55- 

43- 
34- 
26- 



170 

130 
95- 
72- 
55- 

43- 
34- 

26- 
17- 



170 
130 
95 



NSMCE1 



RNF151 



170 
130 
95 
72- 
55 
43 
34 
26 



118: 

95- 
72- 
55- 
43- 
34- 
26- 



B 



IB: FLAG 



IB: HA 



170- 
130- 
95- 
72- 
55- 
43- 
34- 
26- 



130- 
95- 
72- 
55- 
43- 
34- 
26- 



118 

95 

72 

55 

43 

34- 

26- 



170- 
130- 
95- 
72- 
55- 
43- 
34- 



I 



170- 
130- 
95- 
72- 
55- 
43- 
34- 



170- 
130 
95- 
72- 
55- 
43- 

34- 
26 



170 
130 
95- 
72- 
55- 
43- 
34- 
26- 



95- 
72- 
55- 
43- 

34- 
26- 



I 



Figure 7 Ligase activity of putative E3s was validated by in vivo assays. (A) HEK 293FT cells were transfected together with FLAG-E3s and 
HA-Ub. After MG-132 treatment, whole cell lysate was subjected to immunoprecipitation with anti-FLAG beads. Poly-ubiquitination was detected 
by Western blotting with anti-FLAG and anti-HA antibody. (B) Protein mixture was subjected to anti-FLAG beads after heating with 1% SDS. 



et al only contained 343 human E3s based on the HECT 
and RING/U-box domains. Therefore, our mining 
should be more specific than Li et al. and also more sen- 
sitive if the same set of domains were used. Moreover, 
we found that Li et al. mined 34 more E3s based on only 
RING and U-box domains than we did. However, a man- 
ual curation of these 34 genes by us indicated that 26 ac- 
tually did not contain the RING or U-box, 3 were 
pseudogenes, 1 was discontinued gene, 3 were E4s, and 
only 1 typical E3 was left. 

The main goal of the present study was to profile the 
expression of E3s in the mouse testis during spermato- 
genesis. Taking advantage of several multiple tissue 
microarray datasets, we rated the expression levels of 
putative E3s as 5 levels (A, P, HP, MS and SP). The rat- 
ing system was based on a supposed normal distribution 
of gene expression in different tissues and the final rate 
was determined by votes from different datasets. There- 
fore, the rating was highly stringent as was confirmed by 



the RT-PCR results. Almost all predicted testis-specific 
E3s were indeed specifically expressed in the testis. In 
addition, we found that majority of the putative E3s were 
expressed in the mouse testis, and the number of testis- 
specific E3s was much higher than in other tissues such 
as the liver, the lung, and the muscle. The stages-specific 
expression of E3s during mouse spermatogenesis indi- 
cated that most of the testis-specific E3s were expressed 
in spermatocytes and spermatids. As spermatocytes 
undergo meiosis and spermatids go through a number 
of unique processes such as acrosome formation, histone 
replacement, flagellum formation and shedding of most 
of the cytoplasm, it is reasonable that these cells need 
many tissue-specific E3s to catalyze the turnover of 
proteins. 

Conclusions 

We performed a mining of ubiquitin ligases in the 
mouse genome and analyzed their expression profile in 
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the mouse testis systematically. The identification of a 
large number of homologous putative E3s between 
human and mouse suggests that these enzymes are 
highly conserved in the two species. Our analysis reveals 
that 340 putative E3s are expressed in the mouse testis 
and 73 E3s are highly or specifically expressed. Based on 
in vitro and in vivo essays, the ligase activity of several 
novel E3s expressed during the mouse spermatogenesis 
are confirmed. The lists of E3s expressed at different 
levels during spermatogenesis provide a valuable re- 
source from which key factors regulating spermatogen- 
esis could be identified. 

Methods 

mRNA expression evaluation in tissues 

In order to summarize the expression levels of a gene in 
a particular tissue (t) from different datasets (d), the ex- 
pression value of a gene in a tissue (x t ) were converted 
into the z-score (z t = (x t - u)/s where u and s are the 
mean and standard deviation of x t among tissues). As- 
suming the z t values of a gene in different tissues are 
normally distributed, three threshold values (z 2 = 1, 
z 2 = 1.96, z 3 = 2.58) mean that at least 68%, 95%, 99% of 
tissues have the absolute values of their z t scores to be 
less than these value, respectively. For example, if a gene 
in the testis has its z t=testis > 2.58, we can say that the 
z t^testis values of this gene in at least 99.5% of the tissues 
should be less than z t=testis while z t ^ test is values in at 
most 0.5% of the other tissues may be equal to or more 
than z t=testis . The following five expression levels were 
defined accordingly. A for absent of a gene in t if the ori- 
ginal Affymetrix P/A call value is A. P for present if the 
Affymetrix P/A call value is P. HP for highly present 
with Zt > Zi. MS for multi-tissue specific with z t > z 2 be- 
cause a gene with such an expression value may only 
expressed in several tissues. SP for specific with z t >z 3 
for tissue t and Zi < z 2 for all the other tissues (i ^ t). The 
most specific expression level was assigned to a gene if 
it satisfying more than one level. For example, a gene 
considered SP must also satisfy the MS condition, but 
we labeled it as SP. We considered five datasets (includ- 
ing the four mouse microarray datasets and the UniGene 
EST dataset) as five voters. We determined the mRNA 
expression level by a ballot. The expression level that 
was supported by the highest number of votes was 
assigned to the gene in the specified tissue. However, if 
two expression levels were both supported by the high- 
est number of votes, the less specific level was assigned. 
For example, MS was assigned to a gene if both MS and 
SP were supported by 2 votes. 

Animal use 

Male CD1 mice were purchased from Vital River La- 
boratories (Beijing, China) and maintained in the 



Experiment Animal Center, Chinese Academy of 
Sciences. Animal use was approved by the Animal Re- 
search Committee of the Institute of Zoology, Chinese 
Academy of Sciences (Protocol No. 2004-35). Protocols 
for animal sacrifice and tissue harvesting were in accord- 
ance with the NIH Guide for the Care and Use of La- 
boratory Animals. All the mice were euthanized by 
cervical dislocation before tissue harvesting. 

Isolation of different type germ cells 

8 dpp and 17 dpp mice were used for isolating type A 
spermatogonia (SG-A) and pachytene spermatocytes 
(pacSCs), respectively. Meanwhile, adult mice were used 
for isolating round spermatids (rSTs) and elonged sper- 
matids (eSTs). The seminiferous tubules were isolated 
from decapsulated testes and digested into cell suspen- 
sion with collagenase (1 mg/ml, 2 min, 37°C, Gibco) and 
trypsin (0.25%, 5 min, 37°C, Gibco). The purity of four 
types of germ cells all exceeded 90% using the unit grav- 
ity sedimentation procedure in 2-4% BSA medium as 
described before [31-33]. The purity of cells was vali- 
dated based on morphological evaluation and confirmed 
by RT-PCR of reported marker genes expressed in differ- 
ent germ cells. 



Table 4 Primers used in RT-PCR analysis 


Name 


Product length 


Sequence (5'-3') 


Action 


480 bp 


Sense: ATATCGCTGCGCTGGTCGTC 
Antisense: ATATCGCTGCGCTGGTCGTC 


March11 


321 bp 


Sense: ATGAGCGACGAGGGCAAAAAAC 
Antisense: CTACTGAACmCTCAACCAGT 


Nsmcel 


830 bp 


Sense: ATATGGOTCCGGCGTGATC 
Antisense: TOTGCTGCCGGGTCCGTAA 


Pcgf6 


698 bp 


Sense: AAGCCTGCTGCTCCACA 
Antisense: TGACCCACCCACCTCAA 


Rnf126 


924 bp 


Sense: ATATGGCCGAGGCGTCGCCG 
Antisense: TOGAGTOOTGTGGCG 


Rnf19b 


407 bp 


Sense: AGCGGAAGCCCTACAGA 
Antisense: CCACCACCAGCACCA^ 


Rnf151 


519 bp 


Sense: GCCGCCTACAAGTCAAG 
Antisense: GCCTGGGTOAGTGmA 


TrimLI 


1850 bp 


Sense: TGATGTCCAACCATGAGA 
Antisense: CAGATGTGGCAGATGATAAG 


Trim42 


816 bp 


Sense: ACCAATCTGAAGTGCTACTA 
Antisense: CTACAGGTCATOTCAAT 


Trim69 


1517 bp 


Sense: CAAGmCATGGAGGTATC 
Antisense: ATCTGTGGATGTACGATGTG 


Znrf4 


876 bp 


Sense: CATGGCGCGGmGCGTG 
Antisense: TCAGGAGAGCTCGGAAGTCGC1TC 
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RNA isolation and RT-PCR 

Total RNAs were extracted from mouse testes at differ- 
ent developmental stages and other tissues using TRIzol 
solution. 2 ug RNA was reverse -transcribed into cDNA 
by reverse transcriptase (Promega). The primers used for 
PCR are listed in Table 4. PCR reactions were conducted 
following stand protocol Briefly, the reactions were 
started at 94°C for 3 min, and went through 27 cycles 
with denaturing at 94°C for 30 sec, annealing 60°C (or 
other suitable Tms) for 30 sec, and elongation at 72°C 
for 40 sec. Reactions were finished with a final extension 
step at 72°C for 10 min. 

Molecular cloning of putative E3s 

cDNAs of putative mouse E3s were amplified using 
PCRs. Restriction enzyme sites were included in the pri- 
mers. PCR products and vector plasmids were digested 
with selected restriction enzymes, target fragments were 
gel purified and ligated subsequently. Full- or partial 
length cDNAs of putative E3s were sub-cloned into 
pFLAG-CMV4, pEGFP-Nl and pGEX-4Tl to make dif- 
ferent fusion proteins. The primers used for plasmid 
construction were listed in Table 5. All the recombinant 
plasmids were confirmed by DNA sequencing. Cloning 
procedures followed standard protocols. 



Preparation of recombinant proteins 

Putative E3s with GST tag and E2s (UBE2D2, UBE2D3, 
UBE2L3, Ube2k) with His tag were produced in the 
Escherichia coli strain Rossetta (DE3) (TransGen). Bac- 
terial culture, IPTG treatment and lysate sonication were 
conducted following standard protocols. GST and His 
tag fusion proteins were purified by AKTA™ purifier 
9500 using GST-Trap_5ml_FF and His-Trap_5ml_FF 
columns (GE Health care), respectively. Protein purifica- 
tion followed the manufacturers instructions. Purified 
proteins were dialyzed by ultrafiltration (MICROCON, 
Millipore) before in vitro ubiquitination assays were 
performed. 

In vitro ubiquitination assays 

In vitro ubiquitination reactions were set up by adding to 
the Eppendorf tubes with the following reactants: 0.25 \ig 
yeast El (Boston Biochem); 2 \ig various purified recom- 
binant E2s (UBE2D2/UBCH5B, UBE2D3/UBCH5C, 
UBE2L3/ UBCH7, Ube2k/E2-25K); 5 ug putative E3s with 
GST tag; 5 ug MYC-ubiquitin (Boston Biochem) in 
50 mM Tris/HCl (pH 7.4), 5 mM MgCl 2 , 1 mM DTT; 
2 mM ATP (Sigma). Negative control reaction was set up 
with putative E3 being replaced with equal amount of 
glutathione S -transferase. After 30 min incubation at 37°C, 



Table 5 Primers used for recombinant plasmids construction 



Name Vector Sequence (5'-3') 



Nsmcel 


pGEX-4 T1 


Sense: CG GGATCC ATGGOTCCGGCGTGATC 
Antisense: CG GAATO CTAGTGCTGCCGGGTCCGTAA 




pEGFP-N1 


Sense: CG GAATTC ATATGGC1TCCGGCGTGATC 
Antisense: CG GGATCC TOTGCTGCCGGGTCCGTAA 




pFLAG-CMV4 


Sense: CG GAATTC TATGGOTCCGGCGTGATC 
Antisense: CG GGATCC CTAGTGCTGCCGGGTCCGTAA 


Rnf1 26 


pGEX-4 T1 


Sense: CG GGATCC ATGGCCGAGGCGTCGCCG 
Antisense: CG GAATO TCAGGAGTOOTGTGGCG 




pEGFP-N1 


Sense: CG GAATO ATATGGCCGAGGCGTCGCCG 
Antisense: CG GGATCC TOGAGTOOTGTGGCG 




pFLAG-CMV4 


Sense: CG GAATTC TATGGCCGAGGCGTCGCCG 
Antisense: CG GGATCC TCAGGAGTOOTGTGGCG 


Rnf151 


pGEX-4 T1 


Sense: CCG GGATCC TOCTGTGTOTGTCTGCCA 
Antisense: GC CTCGAG COTCATAAGACACmGGC 




pFLAG-CMV4 


Sense: CG GAATTC GATC ATGAGTGGTGG 
Antisense: CG GGATCC TCATAAGACACmGGC 


Znrf4 


pEGFP-N1 


Sense: CCG GAATO CCATGGCGCGGmGCGTG 

Antisense: CCG GGATCC CTGGAG AGCT CGGA AGT CGC TTC 




pFLAG-CMV4 


Sense: CCG GAATO CATGGCGCGGmGCGTG 
Antisense:CCG GGATCC TCAGGAGAGCTCGGAAGTCGC TTC 
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reactions were stopped by treatment at 95°C for 10 min 
in a SDS-PAGE sample buffer containing p-mercap- 
toethanol Then proteins were separated by electrophor- 
esis and subjected to Western blot using anti-ubiquitin 
antibody [34]. 

Cell culture and transfection 

HEK 293FT and CHO cells were cultured in high glu- 
cose DMEM (Gibco) supplemented with 10% fetal bo- 
vine serum (Hyclone) in 5% C0 2 at 37°C 24 hours after 
passage, 3.5 x 10 4 /cm 2 cells were plated and transient 
transfected with Mega Trans 1.0 (ORIGEN) according to 
the manufacturers recommendations. 24 hours after 
transfection, CHO cells were treated with 10 \iM MG- 
132 (Z-Leu-Leu -Leu-al, Sigma) for 6-8 hours before 
fixed in 4% paraformaldehyde for 15 min at room 
temperature. Then the cells were permeabilized with 
0.1% Triton X-100 and washed by PBS. Cell nucleus was 
dyed with DAPI (1 (ig/ml) for 10 minutes. Samples were 
observed with Zeiss LSM710 confocal microscope. All 
the images were quantified using the ZEN software 
(Zeiss). 

In vivo Ubiquitination assays 

HEK 293FT cells transfected with different FLAG-tagged 
putative E3s and HA-tagged PRK-Ub were cultured for 
36 hours. Cells were treated with 10 \iM proteasome in- 
hibitor MG-132 for 6-8 hours. Subsequently, the cells 
were washed twice in ice-cold PBS and suspended in cell 
lysis buffer [50 mM Tris HC1 (pH7.4), 150 mM NaCl, 
1 mM EDTA, 1% NP-40, and 10 \iM protease Inhibitor 
Cocktail (Sigma)]. The cell lysate (3 x 10 7 cells/ml, ap- 
proximately 0.5 ~ 1 mg/ml) was incubated on a shaker 
for 30 minutes at 4°C and centrifuged at 15000 g for 
15 min. 1 ml supernatant was incubated with 30 \A sus- 
pended ANTI-FLAG M2 affinity gel (Sigma) overnight 
at 4°C and then spun at 15,000 g for 15 min [35]. The 
M2 resin was washed 4 times with buffer A [50 mM 
Tris-HCl (pH7.4), 150 mM NaCl], bound proteins were 
eluted using buffer B (including 150 ng/|il 1 x FLAG 
peptide, Sigma). Eluted proteins were separated by elec- 
trophoresis using 4-20% gradient PAGE gels [36,37]. 

In order to disrupt the non-covalent interactions, cells 
were lysed with lysis buffer containing 1% (w/v) sodium 
dodecyl sulfate (SDS) and then heated to 100°C. 10 mM 
sulfhydryl alkylating agent N-ethylmaleimide (NEM, 
Pierce) was added into the cell lysis buffer to prevent 
ubiquitin chain from breaking. The concentrations of 
SDS in the cell lysis were diluted to 0.1% before immu- 
noprecipitated with the ANTI-FLAG M2 beads [38]. 

Western blot analysis 

Western blot analysis followed standard protocol. Primary 
antibodies and dilution factors were as the following: anti- 



FLAG and anti-HA, 1:5000 (mouse, monoclonal antibody, 
Abmart); anti-ubiquitin, 1:1000 (mouse, polyclonal anti- 
body, Santa Cruz). The secondary antibody was anti- 
mouse IgG (goat polyclonal antibody labeled with Dylight 
800, 1:5000, KPL). Fluorescent signal on NC-membrane 
was detected with an Odyssey infrared imaging system 
(LI-COR Biosciences). 
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